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Abstract We discuss three topics related to the neutron 
star (NS) mass spectrum. At first we discuss the possi- 
bility to form low-mass {M ^ 1 Mq) and suggest this 
is possible only due to fragmentation of rapidly rotat- 
ing proto-NSs. Such low-mass NSs should have very high 
spatial velocities which could allow identification. A crit- 
ical assessment of this scenario is given. Secondly, we dis- 
cuss mass growth due to accretion for NSs in close binary 
systems. With the help of numerical population synthe- 
sis calculations we derive the mass spectrum of massive 
(M > 1.8 A/q) NSs. Finally, we discuss the role of the 
mass spectrum in population studies of young cooling 
NSs. We formulate a kind of mass constraint which can 
be helpful, in our opinion, in discussing different com- 
petive models of the thermal evolution of NSs. 
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1 Introduction 



Mass is one of the key parameters for neutron star (NS) 
physics and astrophysics. It can be measured with high 
precision in binary radio pulsar systems. Until very re- 
cently all estimates we re obtained in the very narrow 
region of 1.35-1.45 Mq l|Thorsett fc Chakrabartvlll9 991 
These values lie very close to the Chandrasekhar limit for 
white dwarfs. Thus, M = 1.4 Mq was considered to be 
the standard value of the NS mass. Recently, the range 
widened towards lower masses after the discovery of the 
double pulsar J0737-3039 (Burgav et al. 200 ^1. One of 
the N Ss in this system has M — 1.25 Mq l(Lvne et alJ 
I2fl04l) . Several other examples are known (see Sec. 4). 

The mass range is extended also towards higher masses, 
although these results are less certain. There is one NS 
in a binary radio pulsar system with mass significantly 
higher than the canonical value 1.4 M©. It is the pul- 
sar J0 751-H807 with the mass 2.1 +°1 (95% confidence 
level) (iNice. Solaver fc Stair£ll20d5|] . A recent reanalysis 
of RXTE data suggests that in low-mass X-ray bina- 
ries there are several candidates for neutron stars w ith 
masses between 1.9 and 2.1 Mq ijBarret et al.ll2005l ) as, 
e.g., 4U 1636-536. The small number of massive radio 
pulsars, however, can be a result of selection effects. 

Cooling curves of NSs are strongly dependent on the 
star mass, the mass spectrum is one of the most im- 
portant ingredients of the population synthesis of these 
sources and, unfortunatel y, one of the less known. Still, 
theoretical considerations l|Wooslev. Heger fc Weaveil20n2() 
can shed light on the general properties of the mass spec- 
trum of NSs. This properties have to be taken into ac- 
count when confronting theoretical models of the thermal 
evolution of NSs with observational data. 

In this paper we discuss three issues: the formation 
of low-mass NSs, the mass growth of NSs in binaries, 
and the role of NS masses in the population synthesis of 
young cooling objects. 
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2 Formation of low-mass neutron stars 

The results presented in this section are based on the 
research note (fPoDOv 2004). 

In many models of thermal evolution of compact ob- 
jects (neutron stars - NS, hybrid stars - HyS, strange 
stars, - SS) low-mass sources with 0.8 Ad^ < M < 1 Mq 
^ remain hot for a relatively long time (about few million 
years) ^. During al l that time they remain hotter than 
more massive stars l|Blaschke. Grigorian &: Voskresenskvl 
l2004l and references therein). In that sense they are promis- 
ing candidates to be observed as coolers, and their de- 
tecti on is of great interest for the physics of dense mat 
ter l|Carriere. Horowitz fc Piek arewicz' 



tional waves emission, and after reaching the minimum 
mass 0.1 Mq) it explodes. So, the remaining compact 
object would also have relatively low mass (~ 1 M©) and 
high velocity. 

Objects formed after fragmentation have particular 
predictable properties: high spatial velocity, high surface 
temperature, velocity vector nearly perpendicular to the 
spin axis (since in this mechanism the kick is always 
obtained in the orbital plane which coincides with the 
equator of the initial proto-NS). 

Due to high kicks low-mass compact objects are not 
expected to appear in binaries (at least they should be 
1 '2003") However rarely found in binary systems). To find a low-mass com- 

in m ost of models of NS formation ( Wooslev. HegerTweayp ^ ojjf ^ ^^a^^^h ^""^ ^ ^o*' ^^§1^ 

2nn2 tlFrver fc Kalogeral2nniUTimmes. Wooslev fc Weavei< ^"^^^^ 

199(il) no objects with M < 1-1.2 Mq are formed. It is It is reasonable to expect that mass and kick velocity 
likely that just because masses of stellar cores are always are anti-correlated, as far as a higher mass of the remain- 
heavier than 1.2 Mq even for t he solar metallicity (and ing object corresponds to the lighter exploded compo- 
heavi er for lower metallicity, see lWooslev. Heeer fc Weaveil nent, which means to a wider orbit, and to lower orbital 

velocity of the remaining more massive component. Also 
higher kicks lead to smaller fall-back (Colpi, Wasserman 
2002). 

To reach fragmentation conditions (the dynamical 
instability) it is necessary that the progenitor core is 
rapidly rotating. Rotation of isolated progenitors and its 
influence on properties of newborn NSs was studied in 
several papers (see, for example, Heger et al. 2003 and 
references therein) . To obtain a rapidly rotating compact 
object it is necessary to avoid spin-down influence of the 
magnetic field, so probably compact objects born after 
fragmentation should be low magnetized. It means that 
low-mass neutron stars are not expected to be normal 
radio pulsars. Because of the same reason they are not 
expected to show any kind of magnetar activity. 



Hool. 

In our opinion the only way to form a low-mass object 
from a (relatively) high -mass cor e is fragmentation (see 
however a discussion in IX^l200l where the author dis- 
cusses the formation of low-mass SSs from white dwarfs 
via accretion induced collapse). 

The fragmentation of a rapidly rotating proto-NS 
due to dynamical instabilities as part of a two-stage su- 
pernova (SN) explosion mechanism was suggested by 
iBerezinskii et aL* later the mechanism was de- 

veloned bv llmshennik fc NadezhiiJ (^23) • This mecha- 
nism can explain several particular features of SN explo- 
sions (for example the delay in the neutrino signal from 
SN1987A). In this scenario a compact object inevitably 
obtains a high kick velocit y. Recently, the mechanism w as 
studied in some details bv lColpi fc Wassermanl (|2002tl in 
application to kick phenomenon. Because of high tem- 
perature the exploding (lighter) companion can be sig- 
nificantly more massive than the minimum mass for cold 
stars (^ 0.1 M0), up to ^ 0.7 M©. In that case from 
the initial object of, say, 1.2 - 1.3 solar masses due to 
fragmentation and explosion of the lighter part, we can 
finally obtain a high-velocity NS with a mass of about 
0.8-1 Mq or lower if the mass of the initial object was 
smaller. 

Even lower masses can appear if in the fragmentation 
three bodies are formed. In such a case the lightest or an 
intermediate mass fragment can be dynamically ejected 
from the system (again with significant velocity about 
thousands km s~^). Such ejected compact objects can 
have masses about 0.2-0.5 solar masses. In the remaining 
pair the lighter one can start to accrete onto the second 
companion because of the orbit shrinking due to gravita- 

^ Here and below speaking about compact objects we mean 
the gravitational mass. 

^ Object s with even lower mass^ ~ 0.5 Mr, also a re rela- 
tively hot IBlaschke. Grigorian fc Voskresenskvl l 2OOD . 



We have to note, that the mechan ism of SN explosion 
suggested bv lBerezinskii et all l|l988|) and'imsh enniki |)1992() 
has its internal problems. There is the possibility, that a 
rapidly rotating proto-NS can just loose part of its mass 
in the form of an outflow in the equatorial plane. In that 
case two spiral arms appear, no second (or third) compo- 
nent is formed, and the kick can be relatively small. The 
fraction of lost mass is very small, about 4% (Houser et 
al. 1994), so that the final mass cannot be much lower 
than the initial one. If the fragmentation in the process 
of NS formation never happens in nature, then, in our 
opinion, it is very improbable, that low-mass compact 
objects can exist. The discovery of a high velocity low- 
mass NS, HyS or SS will be a strong argument in favour 
of the mechanism. 

To conclude: the fragmentation of a proto-NS can 
be a unique mechanism of the formation of low-mass 
compact objects, which are expected to have several pe- 
culiar characteristics that can help to distinguish them 
among possible candidates. However, the realization of 
this mechanism in Nature is not very promising. 
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3 Massive neutron stars in binaries 



In this section we present our recent calcu lations of the 
mass growth of NSs in close binary systems l)Popov k, Prokt ~ 
I200fj) . 

Observationally, high masses of NSs are mainly sup- 
ported by data on X-ray binaries, where recently new 
results have been extracted from a reanalysis of RXTE 
data which has suggested that sharp and reproduceable 
changes in QPO properties are related to the innermost 
stable circular orbit. A mass e stimate of 1.9-2.1 M p^ has 
been given for 4U 1636-536 ()Barret et al.ll20nfj) . Esti- 
mates for several systems give also very high mass values 
with rather large uncertainties: 1.8-2.4 Mq for Vela X- 
1^ (Qu aintrell et al. 2003.) 2.4±0.27 Mr., for 4U 1700-37 
Jciark et al. 200 3 see alsolHeineke et al.l2003H van Kerkwii I 




2.6 2,8 
M(NS)/Mo 



l2004|) . IShahbaz et alJ l)2004l) presented observations of a 
low-mass X-ray binary 2S 0921-630/ V395 Car for which 
the l-cr mass range for the compact object is 2.0-4.3 Mq. 

The existence of NSs with M ~ 1.8-2.4 Mq is not 
in contradiction with the present day theory of NS in- 
teriors. There are several models with stiff equation of 
state which allow the existence of NSs wit h masses larger 
than 2 jVfin (see a rev iew and references in lHaense1ll20M 
iLatt imer fc Prakashll20(ifil and in the contribution of J. 
Lattimer in these proceedings). Here we will assume that 
masses of NSs with extreme rotation can reach the max- 
immnvalue_M„iax = 3.45 Mq according to the estimate 
bv lOuved fc 'Butlcr (1999|. If, however, in nature M^ax 
is lower than this value, the part of the obtained mass 
spectrum exceeding it should be attributed to black holes. 



Fig. 1 Mass distribution of NSs. As we are interested only 
in the massive population we do not show the results for 
compact objects with M < 1.8 Mq. The upper mass limit 
corresponds to SkyS with maximum rotation (Ouyed 2004). 
The dashed line represents results for the scenario with zero 
kick. The solid line is the non-zero kick. Left peaks for both 
distributions correspond to NSs with a single episode of ac- 
cretion. Right peaks are formed by NSs which also increased 
their masses via accretion from WDs. Distributions were nor- 
malized to unity, ie. the area below each line is equal to one. 



3.2 Mass spectrum 



As we are interested here in systems with a high mass 
ratio it is necessary to consider three different situations 
after NS formation when the secondary fills its Roche 
3.1 Model lobe: i.) a normal star can fill its Roche lobe without 

common envelope formation; ii.) a normal star can fill 

The m odel we use is discussed in detail inlPopov k ProkhoroiU Roche lobe with common envelope formation; in.) a 
l|2005^ . For our calculations we use the 'Scenario Ma- WD fills its Roche lobe. 



chine" code developed at the Sternberg Astronomical 
Institute. * A description of most of the parameter s of 
the code can be found in Lipunov, Postnov & Prokh orovl 



19961) . The main assumptions of the scenario are: 

— All NSs are born with M = 1.4 A/q. 

— At the common envelope stage a hypercritical accre- 
tion (with M much larger than the Eddington value) 
is possible. 

— During accretion the magnetic field of a NS decays 
to a value that cannot prevent rapid (maximum) ro- 
tation of the NS. 

— The Oppenheimer-Volkoff mass of a rapidly rotating 
NS (the critical mass of a B H formation) is assumed 
to be 3.45 M© according to lOuvedl l)2004|) . 



This range is bas ed on the two estimates given in 
IQuaintrell et aT (lOOTl: 1.88 ± 0.13 and 2.27 ± 0.17 M©. 

http://xray.sai.msu.ru/sciwork/scenario.html and 
/xray.SM.msu.ru~~mystery/articles/review71 



^ http: 
|http://xi 



To fill the Roche lobe a normal secondary star has to 
evolve further than the main sequence stage. During its 
evolution prior to the Roche lobe overfiow the mass of the 
star is nearly constant. A common envelope is not formed 
if the normal star is not significantly heavier than the NS. 
In this regime, mass is not lost from the binary system. 
For more massive secondaries, formation of a common 
envelope is inevitable, mass transfer is unstable. In this 
regime, a significant fraction of the mass fiow is lost from 
the system, so the mass of the NS grows less effectively. 
After the common envelope stage the orbital separation 
becomes smaller, so later on even a degenerate core of 
the secondary - a WD - can fill the Roche lobe. 

In the figure we present our calculations of the mass 
spectrum of massive NSs. About 25% of accreting mas- 
sive NSs have normal stars as secondaries, the rest have 
WD companions. The formation rate of massive NSs was 
found to be 6.7 10~^ yrs~^. This corresponds to ~ 10 000 
of these compact stars in the Galaxy. 
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EqoI. References 
Cooling curves of NS are strongly mass dependent 

(see contributions by D. Page, A. Kaminker and others in 

this volume). Unfortunately, a mass determination with 

high precision is available only for NSs in binary systems. 

Compact objects in X-ray binaries could accrete a sig- 
nificant amount of matter. For some of the radio pulsars 

observed in binaries, accretion also played an important 

role. Without any doubts masses of millisecond pulsars 

do not represent their initial values. However, there is a 

small number of NSs with well determined masses, for 

which it is highly possible that these masses did not 

change significantly since these NSs wer e born (data on 

NS masses can be found, for example, in iLorinie 

and references therein). These are secondary (younger) 

components of double NS systems. According to stan- 
dard evolutionary scenarios these compact objects never 

accreted a significant amount of mass. Their masses lie 

in the narrow range 1.18 - 1.39 Mq. Now there are nine 

double-NS systems with well-estimated masses of the 

secondary components. This set of data is a very good 

evidence in favour of the mass s pectrum used in our 

population synthesis calculations ( Po pov et al.l200(^ . Of 

course, some effects of binary evolution can be impor- 
tant, and so for isolated stars (or stars in very wide bi- 
naries) the situation can be slightly different. However, 

with these observational estimates of initial masses of 

NSs we feel more confident using the spectrum with a 

smaU number of NSs with M ^ 1.4 - 1.5 Mq. 

Brighter sources are easier to discover. So, among 

known cooling NSs the fraction of NSs with typical masses, 

i.e. in the range 1.1 M© ~ M ~ 1.5 M©, should be even 

higher than in the original mass spectrum. So, we have 

the impression that it is necessary to try to explain even 

cold (may be with an exception of 1-2 coldest) sources 

with M ~ 1.4 - 1.5 Mq. Especially, the Magnificent seven 
and other young close-by compact objects should be ex- 
plained as most typical representatives of the whole NS 
population. We want to underline that, even being se- 
lected by their observability in soft X-rays, these sources 
form one of the most uniform samples of young isolated 
NSs. In this sense, the situation where a significant num- 
ber of sources are explained by cooling curves corre- 
sponding to 1.5 Af© ^ M ^ 1.7 Mq should be consid- 
ered as a disadvantage of the model. Particularly, Vela, 
Geminga and RX J1856-3754 should not be explained as 
massive NSs. 

All the above gives us the opportunity to formulate 
the conjecture of a mass spectrum constraint: data points 
should be explained mostly by NSs with typical masses. 
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